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Abstract. The precursor ESA mission LISA-Pathfinder, to be flown in 2013, aims at 
demonstrating the feasibility of the free-fall, necessary for LISA, the upcoming space-born 
gravitational wave observatory. LISA Technology Package (LTP) is planned to carry out a 
number of experiments, whose main targets are to identify and measure the disturbances on each 
test-mass, in order to reach an unprecedented low-level residual force noise. To fulfill this plan, 
it is then necessary to correctly design, set-up and optimize the experiments to be performed 
on-flight and do a full system parameter estimation. Here we describe the progress on the non- 
linear analysis using the methods developed in the framework of the LTPDA Toolbox, an object- 
oriented MATLAB Data Analysis environment: the effort is to identify the critical parameters 
and remove the degeneracy by properly combining the results of different experiments coming 
from a closed-loop system like LTP. 



1. Introduction 

LISA [1] will be the first gravitational wave detector in space. Together with the development 
of the flight hardware, like the Gravitational Reference Sensors (GRSs) and the interferometers, 
data analysis plays a crucial role. In order to prove the feasibility of the free-fall, a scientific 
mission, LISA-Pathfinder (LPF) [2,3], has been developed. To date, the mission is in the final 
stage before the launch: all the subsystems are being fully tested and mostly assembled [4]; at 
the same time, data analysis is carrying on [5-7]. 

On board LPF, the LTP experiment [8] will verify the free-fall with a unprecedented low- level 
in force disturbances. In order to achieve this, a measurement scheme has been proposed: two 
Au-Pt 2-kg Test-Masses (TMs) separated by ~30 cm reproduce one Doppler-link among the 
six within LISA. A laser interferometer keeps track of the relative displacements and rotations 
between one TM and the spacecraft and between the two TMs. Also, two GRSs provide the 
locking of the TMs along the non-optically sensed axes. The aim is to measure the residual 



forces or any systematics on the first TM, subtract them from the data and improve the quality 
of the free-fall. A key technique consists on the drag-free: the TM is let in free-fall and the 
spacecraft is forced to follow the TM itself, based on a suitable readout. At the same time, the 
second TM is forced to follow the first, in an analogous scheme. The requested forces/torques 
on the spacecraft and the second TM are then calculated by the Drag- Free Attitude and Control 
System (DFACS). 

In order to measure the force disturbances on the TMs, it is necessary to implement a set 
of identification experiments, whose main targets are to extract information on the system 
critical parameters. Each experiment is specifically designed to excite the system and hence its 
parameters at proper frequencies. The modeled system output, compared to the data, can give 
estimates for the parameters, but unfortunately some of them will be degenerate. However, by 
using all the system outputs at the same time, the information is maximized and the parameter 
degeneracy can be removed: in the end the system is fully identified. 

In this paper it is shown how to perform a non-linear time-domain parameter estimation of 
mission-like simulator data during an operational exercise. In Section 2 we describe the physical 
model; in Section 3 we introduce the analysis technique that will be applied subsequently to 
mission-like simulator data in Section 4. Final comments are discussed Section 5. 

2. LTP as a closed-loop physical system 

As introduced previously, the main LTP science mode consists on one TM in drag-free along the 
laser-sensitive axis (x-axis) and suspended by electrostatic actuators along the other degrees of 
freedom. The DFACS, based on the interferometer read-out, commands the spacecraft to follow 
the first TM and the second TM to follow the first. 

The 1-D dynamics for the two TMs and the spacecraft (SC) can be easily written in frequency 
domain as 

mi s 2 x\ = F\ — m\ T x ^ 2 (x2 — xi) — mi u>fx\ — mi s 2 xsc > 
m 2 s 2 x 2 = F 2 + m 2 T Xy i (x 2 — xi) - m 2 oj 2 x 2 - m 2 s 2 xsc 

+ m 2 A sus H sus oi 2 , 
msc s 2 x sc = F S c + mi ivfxi + m 2 u\x 2 - F 1 - F 2 

- m 2 A sns H sus oi 2 + m S c A df H di oi , 

where s is the Laplace frequency and xi and x 2 are the two TMs' coordinates with respect 
to a local non-inertial reference frame, fixed to the spacecraft. In Eq. (1), Fi, F 2 and F$c are all 
the external forces on both TMs and spacecraft; mi ufxi and m 2 co 2 x 2 are the residual coupling 
modeled as oscillator-like forces, being oj\ and uj 2 named the parasitic stiffnesses; T Xt i and T x ^ 2 
are the mutual gravity gradients between the TMs. The first loop term is msc Aif -f^df °i which 
is the DFACS commanded force actuated on the spacecraft by the micro-propulsion thrusters, 
based on oi, the interferometer reading of x\. The second loop term is m 2 A SVLS H sns oi 2 , which 
is the DFACS commanded force actuated on the second TM by the electrostatic suspensions, 
based on 012, the interferometer reading of x 2 — x\. Moreover, the loop terms are split into the 
DFACS control laws, and H sus , and the actual application of the commanded forces, A d f 
and ^4 S us 

Since the optically measured coordinates are xi and X12 = x 2 — xi, whilst xsc is n °t directly 
detected, it is natural to eliminate the latter via elementary algebra and define the new physical 
variable x = ). Provided that: 

01 : the interferometer read-out of the displacement Xi of the spacecraft to the first TM; 
012 : the interferometer read-out of the displacement xi 2 of the second TM to the first TM; 



then x converts to the interferometer read-outs through o = S • x + o n , where o = ("12 )> °n 
are the interferometer read-out noises and S is the interferometer conversion matrix containing 
possible cross-coupling terms as off-diagonal terms. Hence, the first two equations of motion in 
Eq. (1) can be rearranged in a single matrix equation and, by solving with respect to o, it is 
possible to demonstrate that o = o s + o n , which contains a deterministic part o s and a noise 
contribution o n from the read-out and force noises. The deterministic part is essentially made of 
injected stimuli to the system both in the controller setpoint (interferometer guidances) and as 
direct applied forces. In this paper we refer to the analysis of only the first case. Consequently, 
if Oj are the inputs to the controller setpoint, the deterministic part is given by o s = H mc Q • Oj, 
where H m( ji is the transfer function matrix model from the input guidances to the interferometer 
outputs. 

3. Parameter estimation technique 

Once derived the deterministic part o s , it is logical to assert that exciting the system in one 
channel and analyzing both outputs, the system can be fully characterized and information on 
the critical parameters can be extracted. Finally, this knowledge takes to the exact computation 
of the residual force noise on each TM [9] , which is a key target of the mission. 

Recall now what is assumed to be known: firstly, the interferometer transfer function matrix 
obtained from the equation of motion, say H m( ji, which depends on p, the set of all critical 
parameters (see Section 4.1); in the end, the interferometers injections Oj. Therefore, the 
modeled outputs to be compared to the interferometer read-outs o^ata hi time domain are 

o md i(i, p) = St 1 [ H mdl(^, P) • Oi(uj)] , (2) 

where Jj" 1 is the inverse Fourier operator and oj the Fourier frequency. The comparison between 
data and model is performed on whitened time-series in a least-square sense. The whitening 
technique is part of a more general framework, the colored noise generation [10]. 
As a brief summary, the estimation: 

(i) assumes gaussianity and stationarity; 

(ii) treats whitened data and models; 

(iii) works in time domain; 

(iv) is non-linear since H mc ji in Eq. (2) is so with respect to each parameter; 

(v) is a multi-experiment/multi-channel analysis. 

4. Analysis of operational exercise data 

The aim of this section is to apply the concepts described in Section 2 and 3 to simulated data. 
The LISA-Pathfinder STOC Simulator (LSS) is a software written for the ESA-STOC center by 
the aerospace industry to fulfill the requirements of very realistic mission-like data and test all 
the methods developed by the analysts. Recently, it has been proposed to simulate a number 
of exercises and practice on these data sets. Here the system identification of an operational 
exercise has been performed in a very mission-like manner. The final goal is to give an estimate 
of the best-fit parameters and their uncertainties, in a similar way it will be done during the 
mission. 

4-1- System parameters 

Here is a description of the system critical parameters being considered in the analysis: 

Htf : drag-free gain; 

H sus : electrostatic suspension gain; 



#21 : interferometer sensing matrix coefficient, named cross- coupling or cross-talk, from on and 
012; 

ujf : stiffness of the first TM; 
Cl>2 '• stiffness of the second TM; 

A^ioop, 1 : time-delay in the control loop for the first read-out; 
A^ioop, 12 : time-delay in the control loop for the second read-out; 
n : characteristic time for the thruster response; 
T2 : characteristic time for the electrostatic suspension response. 

4-2. Stating the experiments 

When defining the type of identification experiments to be performed during the mission, there 
is a general freedom in selecting the injections to be applied to the system. However, the system 
is resolved and the information maximized if the proper frequencies are input: a series of sine 
waves makes the role of the actual stimulus to the system being studied. 

Two general classes of experiments are described: the unmatched- and matched- stiffness. In 
the second case, the control mode is set up so that the stiffnesses are artificially matched to the 
same value. In this paper we focus only to the unmatched-stiffness case and, in order to estimate 
the physical parameters, two types of identification experiments can be thought. Therefore we 
refer to: 

Exp. 1 when injecting 0^1 into the interferometer guidance for 0\; 
Exp. 2 when injecting 0^12 into the interferometer guidance for oyi- 

For the sake of clarity, in Figure 1 the interferometer outputs and guidances are shown for 
both experiments. The inputs that define the two injection experiments are made of a series 
of sine waves of pre-determined frequencies selected to excite specifical elements of the transfer 
function matrix. 
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(a) (b) 
Figure 1. Interferometer outputs 01, 012 and guidances for (a) Exp. 1 (0^,1) and (b) Exp. 2 (0^12) 



4-3. Estimation 

The parameter estimation is performed by following the ideas in Section 3 and the modeled 
outputs are computed through Eq. (2). The whitening filters are computed on long noise data 
stretches before the injection and applied both to the modeled outputs at each iteration step and 
to data. The estimation is focused on a joint analysis of all four read-outs coming from both 
experiments: this maximizes the information and allows the identification of all parameters. 
For a detailed description on non-linear estimation methods and their issues, see for example 



[13]. The algorithm has been implemented in the LTPDA Toolbox [11] under the MATLAB 
programming interface [12]. 

Final results are contained in Table 1. As it is clear, except for H SVLS and S21, whose best-fit is 
compatible (within the error) to the respective initial guess, all the others are moved well apart. 
The errors are computed in the linear approximation, hence only analytical first derivatives 
of the modeled output are taken into account. The best-fit values for this data set are also 



Table 1. Best-fit parameters at x 2 jv = 1.2 and v — 21964. 



Parameter 


Guess 


Best-fit 




1-CJ 


error 


Hdf 


1 


1.07988 




6.5 


10~ 4 
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0.999998 




3.9 


10- 5 


S21 
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-6 


1.2 
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10~ 9 




-1.4- 10" 6 


-2.0356 • 
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2.5 


10 -9 


Atloop, 1 [s] 


0.3 


0.19875 




4.8 


10" 4 


Atloop, 12 [s] 


0.3 


0.2048 




8.0 


10~ 3 


n M 


0.1 


0.4109 




1.9 


10~ 3 


r 2 [s] 


0.01 


0.1984 




2.8 


10~ 3 



confirmed by other independent estimation methods developed within the LTPDA team, so we 
are confident on the numbers proposed in Table 1. 

Figure 2 shows the plots of the best-fit model compared to the respective read-out, both 
whitened, (panels (a) and (b)) for Exp. 1 and Exp. 2, respectively. The agreement is also 
confirmed by the fact that the final residuals are compatible with a zero-mean gaussian white 
noise with unitary standard deviation. In particular, for Exp. 1 the agreement between data and 
model is acceptable expect for the first read-out (see panel (a) in 2) where some systematics still 
remain, but of the order of 1%. Instead, for Exp. 2 it is clear that the first read-out (see panel 
(c) in 2) contains no signal: indeed there is no cross-coupling S12 between the two read-outs: 
that is why we have not taken into account in the analysis. 

As a final consideration, it is important to point out that the approach to parameter 
estimation in this paper is clearly experimental: no one exactly knows which are the real 
parameters. Anyway, some practical considerations and the agreement with other methods 
suggest that in principle the method is capable to extract the system critical parameters, among 
those, the most important are surely the two stiffnesses and the interferometer cross-coupling. 

5. Final remarks 

The main scientific pay load onboard LPF, the LTP experiment, will demonstrate that is possible 
to keep a TM in free-fall with an unprecedented low level in force disturbances. After the 
completion of the LPF mission, this expertise will be transferred to LISA, the gravitational 
wave observatory in space. The LTPDA team is working in developing all the data analysis 
procedures to characterize and study the system. For it, parameter estimation has a key role 
and in this paper a method has been proposed to be used during the mission. 

The proposed method is a non-linear, time-domain, least-square minimization of a joint 
X 2 that embraces the parameter information from all experiments and read-outs. The overall 
complication is twofold. From one side, it is easy to write the dynamics in frequency domain but 
its impossible to derive analytical models in time domain. From the other side, the LTP noise 
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Figure 2. Whitened best-fit model and data for: (a) Exp. 1 first read-out, (b) Exp. 1 second read-out, (c) 
Exp. 2 first read-out, (d) Exp. 2 second read-out. 



shape is not flat for both read-outs and the whitening filters are derived to properly whiten data 
and models. In the end, the whitening allows us to treat data in least-square sense and compute 
both the best-fit parameters and their covariance matrix. 

The STOC simulator has given the opportunity to use and test the method in the same 
way as it could be done during the LPF mission. The non-linear estimation technique has 
shown appreciable results evident both in the plots of the best-fit curves and the analysis of the 
residuals, but also in agreement with other independent methods. 

Further investigations will surely come up on new STOC data with increasing realism, until 
the mission will start in the next years. For the moment, the time-domain non- linear method 
seems to be one of the most promising parameter estimation solution to be considered for the 
mission. 
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